The Discovery of a Population of Pre-Main Sequence Stars in 
N66/NGC346 from Deep HST/ACS Images^ 



O 

o 

(N 
D 



> 

00 

(N 
(N 
O 

o 

:^ 

o 
U 



X 



A. Nota^'^, M. Sirianni^'f E. Sabbi^ M. Tosif M. Clampinf J. Gallagherf M. Meixner^, S. 
Oey[ A. Pasquali? L. J. Smith? R. Walterbos}°J. Mack^ 



ABSTRACT 

We report the discovery of a rich population of low mass stars in the young, massive star 
forming region N66/NGC346 in the Small Magellanic Cloud, from deep V, I and Ha images taken 
with the HST/ACS. These stars have likely formed together with the NGC346 cluster, ~ 3-5 Myr 
ago. Their magnitude and colors are those of pre-main sequence stars in the mass range 0.6-3 
Mq, mostly concentrated in the main cluster, but with secondary subclusters spread over a region 
across ~ 45 pc. These subclusters appear to be spatially coincident with previously known knots 
of molecular gas identified in ground based and ISO observations. We show that N66/NGC346 is 
a complex region, being shaped by its massive stars, and the observations presented here represent 
a key step towards the understanding of how star formation occurred and has progressed in this 
low metallicity environment. 
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N66 is the most massive and active star form- 
ing region in the Small Magellanic Cloud (SMC) 
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and the star cluster NGC346 is its ionizing core. 
It contains more than half of the known O stars in 
this galaxy, and associated CO clouds (Rubio et 
al. 2000). The upper end of its stellar mass distri- 
bution has been well studied in the recent years: 
very early 0-type stars (03) and high mass (45- 
100 Mq) stars were found suggesting a slope of the 
Initial Mass Function (IMF) in NGC346 very simi- 
lar to the slope for massive stars in LMC and solar 
neighborhood (F = d(log ^) /rf(log M) = -1.9), and 
an age of ~ 3 Myr (see Bouret et al. 2003, and 
reference there in). Massey, Parker & Garmany 
(1989) also suggested presence of sequential star 
formation in NGC346, with star formation begin- 
ning at the southwest side of the association, and 
spreading towards the central cluster. 

While a wealth of information is available on 
the upper end of the stellar mass distribution in 
N66/NGC346, httlc is known of the low end of the 
distribution. The SMC is characterized by a dust- 
to-gass mass ratio that is 30 times lower than in 
the Galaxy (Stanimirovic et al. 2000), and a metal- 
licity about 1/5 than that of the Galaxy. May 
these differences have had an impact on the clus- 
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ter formation mechanism and their distribution of 
stellar masses? To answer this question we began 
a comprehensive investigation of star formation in 
the SMC, from the ground and with HST. This 
paper presents the first results on N66/NGC346. 

2. Observations 

Wc used the HST/ ACS Wide Field Channel 
(WF) to map NGC346 in the filters F555W (~V), 
F814W (~I), and F658N (~ Ha). The observa- 
tions were carried out between July 13 and 17, 
2004, as part of Cycle 13 proposal 10248 (PI. A. 
Nota). In the broad band filters, nine deep images 
(200" X 200" each) were obtained in each filter, 
partially overlapping to get optimal spatial cover- 
age of the NGC346 region. An additional set of 
short exposures was taken in the same instrumen- 
tal configuration to perform accurate photometry 
of the brightest stars. A single image of NGC346 
was taken in the light of Ha, with an exposure 
time of 1542s sec. 

In Plate 1, we show the fully reduced, multidriz- 
zled and combined V-I color composite image of 
NGC 346. The image covers an area of ~ 5' x 5' 
which corresponds to ~ 87 x 87 pc, at a distance 
to the SMC of 59.7 Kpc (Dolphin et al. 2001). The 
exposure time in the central region is ~ 4100 sec 
in each filter. In Plate 2 we show the fully reduced 
Ha image. 

All images were reduced rising the standard 
STScI ACS pipeline CALACS (HST Data Hand- 
book for ACS). 

Photometry was performed on the drizzled 
images by aperture and PSF fitting, using the 
DAOPHOT package within IRAF. Stars were au- 
tomatically detected with the routine DAOFIND 
in all the frames, with a detection threshold set 
at 4iT above the local background level. Features 
misinterpreted for stars, such as noise spikes, were 
rejected. A spatially variable PSF was computed 
for both V and I images, using~ 140 isolated stars 
in different positions of the images. No CTE cor- 
rection was applied. However, preliminary CTE 
determinations indicate that such correction is 
small, and can therefore be neglected in this first 
paper without changing the overall conclusions. 
Our final photometry is calibrated in the Vega- 
mag system, with the zcropoints adopted from 
Sirianni et al. (2005) of 25.724 and 25.501 for the 



V and I filters respectively. 

3. The morphology of the NGC346 region 

The deep HST/ ACS images show an unprece- 
dented level of detail in the NGC346/N66 com- 
plex. The NGC346 cluster at the center is resolved 
into at least three subclusters, and a myriad of 
small compact clusters is visible throughout the 
region. Some of these clusters appear to be still 
embedded in dust and nebulosity, and they might 
be sites of recent or ongoing star formation. The 
diffuse nebulosity is due to the presence of H/3 and 
[OIII] 4959, 5007 in the V filter, which indicate 
the bright underlying HII region. This extends 
to the N with arches and ragged filaments, but is 
delineated to the S-SW by a sharp fan-like ridge, 
which encompasses the entire cluster. This ridge 
is especially evident in the Ha image of the re- 
gion (Plate 2). Extensive dust lanes are well visi- 
ble along this ridge. Dust patches arc; also visible 
throughout the region, most likely the remnant of 
the original molecular cloud (Mizuno et al. 2001). 

4. The Color-Magnitude Diagram 

The V and I photometric data with photometric 
error smaller than 0.1 mag was used to construct 
the V, V-I Color-Magnitude Diagram (CMD) of 
the entire NGC346 region. The CMD (small black 
dots in Fig. 1) contains ~ 80000 stars. Thanks 
to the high photometric accuracy obtained with 
the ACS, these stars produce a very clean CMD, 
with a well defined narrow main sequence (MS), 
down to V 26. The photometric errors in magni- 
tude and color are indicated directly in the CMD, 
as side bars. A pronounced old main sequence 
turnoff is clearly visible, at V ~ 22, and so is the 
red clump, at V ~ 19.5. The CMD morphology 
also reveals a young population, indicated by the 
bright, blue MS, extending from the MS of the 
old SMC stellar population. Below V ~ 21 mag a 
group of very red faint stars is detected to the right 
of the MS. Plausibly, this split MS could be due 
to differential reddening. However, the absence of 
such a split in the upper part of the MS led us 
to believe that these are stars that are less mas- 
sive than ~3 Mq (down to 0.6 M0), which formed 
with the rest of the central cluster but have not 
reached the MS yet (pre-MS stars). 

Analyzing separately the stellar content of the 
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individual clusters wc resolve in the ACS images, 
we can better characterize the local stellar popu- 
lation. The orange dots in Figure 1 highlight the 
observational CMD for the core of NGC346. They 
indicate stars within a radius of ~ 8" from the 
cluster center. This is superimposed to the CMD 
of the entire region (small black dots). 

We superimposed to these CMDs isochrones 
for pre-MS stars from Siess, Dufour & Forestini 
(2000), and for MS and post-MS stars from Lcjc- 
une & Schaerer (2001), and Pietrinferni et al. 
(2004), to define the age and the nature of the 
population observed in the central cluster. We 
find presence of at least two different populations: 

1. a composite SMC field population covering 
a wide range in age, but dominated by the 4 
Gyr component, 

2. a very young population of age ~3-5 Myr. 
Stars less massive than ~3 Mq (down to 0.6 
Mq) are still in the pre-MS phase, approach- 
ing the MS. 

An intermediate age population (~150 Myr) 
might also be present. 

5. The candidate pre-MS stars 

The pre-MS stars are, as expected, faint (V < 
21) and red (1.5<V-I <2.2). The CMD isochrone 
fit indicates they likely formed with the rest of 
the cluster (~ 3-5 Myr ago), but have not reached 
the MS yet. They are in the mass range 0.6 
- 3 M0. Their spatial distribution is shown in 
Figure 2. They appear concentrated at the cen- 
ter of NGC346, mostly in the central subclusters, 
but also in the small compact clusters that are 
located on a direction which is perpendicular to 
the NGC346 main body, extending to the NE. 
Their spatial location appears coincident with the 
clumps of neutral material identified by Rubio et 
al. (2000). The spatial density of the pre-MS stars 
appear to decrease as we move from the center to- 
wards the S-SW ridge, where we observe the dust 
lanes. As also suggested by Rubio et al. (2000), it 
is likely that the region, and especially the ridge, 
may be a site of more recent or on going star for- 
mation, and the forming stars are still embedded 
in dust and too faint to be visually detected in the 
ACS images. 



6. Star formation in NGC346 

The region around N66/NGC346 is a remark- 
able region that has displayed a high level of activ- 
ity in the recent past. To the NE of N66/NGC346 
there is a known SN remnant, SNR0057-7226, first 
identified as a bright source in Einstein X-ray 
observations (Inoue, Koyama, & Tanaka 1983). 
From the expansion velocity derived from FUSE 
observations of ~150 km (Danforth et al. 
2003), we can derive a dynamical timescale of 2 
X 10^ yr and set a lower limit to the epoch of 
the SN explosion. It is clear that SNR0057-7226 
has been interacting with the surrounding region, 
and the O VI CIII and X-ray emission observed 
(Danforth et al. 2003) arises from the shock in- 
teraction. Two concentric bubbles are seen to the 
NE in the Ha region shown in Plate 2: we believe 
they are most likely shaped by the powerful winds 
of HD5980. 

Observations of Rubio et al. (2000) and Con- 

tursi et al. (2000) have painted a detailed pic- 
ture of the molecular gas and dust component in 
N66/NGG346. Their observations indicate that 
the molecular gas has mostly been dissociated by 
the UV radiation of the hottest stars: what has 
not been dissociated is detected in small, com- 
pact clumps, which appear to be spatially con- 
cident with the regions where we observe the pre- 
main sequence stars (Rubio et al. 2000). Specif- 
ically, their C0(2-l) SEST observations indicate 
clumps of emission in the NGC346 central cluster, 
and, even stronger emission along the perpendic- 
ular direction to the NGC346 main body, to the 
NE. These clumps are exactly where we find the 
highest concentration of pre-MS stars. There ap- 
pear not to be much molecular gas or dust outside 
the NGC346 main body, casting the doubt that 
the SW ridge might not be at the boundary with 
a molecular cloud as initially suspected (Rubio et 
al. 2000). 

We have shown that N66/NGC346 is a com- 
plex region that has been shaped by the evolution 
of the most massive stars, some of which have al- 
ready exploded as supernova and are compressing 
the central region from different directions. The 
outstanding question will be to understand how 
star formation has been triggered and has been 
progressing in this region. Further necessary work 
on this region includes spectroscopic characteriza- 
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tion of the prc-main sequence stars, and spectro- 
scopic observations of the gas kinematics. On the 
theoretical side, NGC346 can easily become the 
ideal laboratory to test current theories for star 
formation at low masses. 
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This 2-column preprint was prepared with the AAS lAT^jX 




Fig. 1. — The Color Magnitude Diagram of the 
NGC346 central cluster (large orange dots) super- 
imposed to the CMD of the entire region (small 
dots). Two isochrones (5 Myr and 4.2 Gyr) are 
overplotted for reference. 5 




R.A. 



Fig. 2. — The spatial distribution of the pre-MS 
stars (black dots) shown superimposed to the 7fim 
emission peaks from IS0/LW2 (from Fig. 8 of 
Rubio et al. 2000) 
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Fig. 3. — Plate- 1 VI color composite of three 
partially overlapping HST/ACS WF images of 
NGC346, covering a total area of 5' x 5'. North 
is up and East to the left. The yellow circle 
(radius=8") marks the position on NGC346 
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Fig. 4.— Plate-2 Ha HST/ACS WF image of the 
central NGC346 region, covering an area of 200" 
X 200 ". North is up and East to the left. 
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